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Abstract 

Discrete Mechanics based on finite element methods is presented in this paper. 
We also explore the relationship between this discrete mechanics and Veselov dis- 
crete mechanics. High order discretizations are constructed in terms of high order 
interpolations. 

Keywords. Discrete mechanics, Finite element methods, Symplectic integra- 
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1 Introduction 

We begin by recalling the variational principle for Lagrange mechanics. Suppose Q 
denotes the configuration space with coordinates q % and TQ the tangent bundle with 
coordinates {q l ,q l ) 1 i = 1,2, ••• ,n. In the following, we will work with the extended 
configuration space R x Q with coordinates (i,<7 4 ) and the extended tangent bundle 
R x TQ with coordinates (t,q l ,q l ). Here t denotes the time 0. 

Consider a Lagrangian L : R x TQ — > R. The corresponding action functional is 
denned by 

S((t,q l (t))) = f L{t,q*{t),q\t))dt, (1.1) 

J a 

where q l (t) £ C 2 ([a, &]), which is the set of all C 2 curves in Q denned on [a, b). 

The variational principle seeks the curves q l {t) in C 2 ([a, &]), for which the action 
functional S is stationary under variations of q l (t) with fixed endpoints. 

Consider a vector field on R x Q 

^ = £(W)^+£^WW (i.2) 
f=i H 

Let F e be the flow of V. For (t, q>) £ R x Q, we have F e (t, q>) = (t, q l ). 

t = f(e,t,qi), (1.3) 
<f=.9'(e,i,^), (1.4) 
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where 



d_ 

de 
d 
1c 



f(e,t,<?) =at,Q j ), 



e=0 



(1.5) 



i {e,t,<t) = 4> i {t,(?). 



The transformation (1.31.4) transforms a curve q J (t) in Q into another curve q l (e,t) in 
Q determined by 

t=f(e,t,qi(t)), 
q i =g i (e,t,<?(t)). 

Before calculating the variation of S, we consider the first order prolongation of V 



£3 £5 O 

prV = £(i, <f)|- + E J7 + E ai (*> «^); 



<9i f-^ ' dq 

where pr 1 !/ denote the first order prolongation of V and 

a% qi.qi) = D t <f>\t, q*) - q'D^t, 

where D t denotes the total derivative jl4j]. 

Now we calculate the variation of S directly 



dq* 



(1.6) 



(1.7) 



SS 



de 



c=0 



d_ 

de 



e=0 •/« 

br/ dL d 



L[t 1 q\e,t) 1 j i q\e,t)) dt 



<)L 



dt dt \dq 



ax d <9L 
rft dq 1 



dt 



L 



dL 
dq 1 



(1.8) 



If £(a,g J (a)) = £_(b,q 3 (b)) = and <^'(a,g J '(a)) = ^(b,q j (b)) = 0, the requirement of 
<55 = yields the energy evolution equation 



<9L 



<9L 



^ dt \dq iq 
from £ and the Euler-Lagrange equation 

dL d dL 



L 







5g l di <9g ? 



= 



(1.9) 



(1.10) 



from 

If L does not depend on t explicitly, i.e., L is conservative and ^ = 0, then ( [Tg| ) 
becomes the energy conservation law 



d f dL 



dt \dq 



q l -L\=0 



(1.11) 



1 Note that Hamilton's variational principle holds in fact for the vertical variation, while the horizontal 
variation gives rise to the conservation relation between the Euler-Lagrange equation and energy con- 
servation, see, for example, [H. Here, in order to transfer to discrete case more directly, this variational 
requirement is employed. 
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If wc drop the requirement 



0W(a)) = 0WW) = O, 
we can define from the second term in ( |1.8| ) the extended Lagrangian one form 

We define the extended Lagrangian two form to be 

10 L = d9 L . 



Suppose g l (t,v q ) is a solution of (L1C) depending on the initial condition v q S TQ. 
Restricting q l (e,t) to the solution space of ( 1.10 ) and using the same method in Q, we 



can prove 

(pv 1 g)*w L =w L , (1.13) 

where pr 1 g l (s,v q )) — (s, g i (s,v q ), ■4-g i (s,v q )) denotes the first order prolongation of 
g l {s,v q ). 

By defining a discrete Lagrangian with fixed steps and constructing a discrete action, 
Veselov develops a kind of discrete mechanics [lfj. See also [jlj| [l?]]. The resulting 
discrete Euler-Lagrange equation (variational integrator) preserves a discrete Lagrange 



two form. However the discrete versions of ( [L9| ) and (1.11) are not reflected in Veselov 
discrete mechanics. Using variable steps and a defined discrete energy, Kane, Marsden 
and Ortiz obtain symplectic-energy integrators . Chen, Guo and Wu develop a discrete 
total variation and prove the symplectic-energy integrators preserve a discrete extended 
Lagrange two form || . Guo and Wu also present the difference variational approach with 
variable step-lengths 

The connection between Veselov discrete mechanics and finite element methods was 
first suggested in jl^] . Symplectic and multisymplectic structures in simple finite element 
methods are explored in Q. As is known, it is difficult to obtain high order discretizations 
in Veselov discrete mechanics. We find that the finite element methods provide a natural 
way to obtain a kind of discrete Lagrangian, including those of high order. The purpose of 
this paper is to develop such kind of discrete mechanics based on finite element methods. 

This paper is organized as follows. In the next section, we will present the discrete 
mechanics based on finite element methods. Section 3 is devoted to discussing the role 
of time steps. We finish this paper by drawing some conclusions in Section 4. 



2 Discrete mechanics based on finite element methods 

2.1 Generalized Veselov discrete mechanics 

We first present Veselov discrete mechanics in a more general formulation that allows 
for variable time steps. In this generalized Veselov discrete mechanics, M x M is used 
for the discrete version of R x TQ. Here M = R x Q. A point (t k , q k ; tk+i, <Zfc+i) S 
M x M corresponds to a tangent vector ^"^Zj?* . Here for brevity of notations, we 
restrict ourselves to one-dimensional configuration space Q. A discrete Lagrangian is 
defined to be L : M x M — > R and the corresponding action to be 

AT-l 

S = h(t k ,q k , t k+1 ,q k+1 )(t k+1 - t k ). (2.1) 
fc=o 
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The discrete variational principle is to extremize S for variations of q k holding the end- 
points go and yjv fixed. In this discrete variational principle, t k play the role of parame- 
ters (They can also play the role of variables, see next section). This discrete variational 
principle determines a discrete flow $ : M x M — > M x M by 



$(tfc-l) <7fc-l)*fc) Qk) = (tk,qk,tk+i,qk+i)- 



(2.2) 



Here qk+i is found from the following discrete Euler-Lagrange equation (variational in- 
tegrator) 

(ife+i - tk)D 2 h(t k ,q k ,t k+1 ,q k+1 ) + (t k - t k _ 1 )D 4 h(t k _ 1 ,q k _ 1 ,t k ,q k ) = (2.3) 

for all k e {1, 2, ■ • • , N — 1}. Here Di denotes the partial derivative of L with respect 
to the ith argument. Notice that t k are parameters here. The reason we still call them 
arguments is that we want to keep the notations consistent since in next section they 
will play the role of variables. 

Now we prove that the discrete flow $ preserves a discrete version of the Lagrange 
two form 



del, 



oq 



(2.4) 



We do this directly from the variational point of view, consistent with the continuous 
case Q. 

As in continuous case, we calculate d§ for variations with varied endpoints. 



■ • • , <?jv) ■ {5q ,5qi, ■ ■ ■ , Sq N ) 
§(go + tSqo, qi + eSqi, ■ ■ ■ ,qN+ e<5gjv) 



dS(q ,qi, 

d_ 
~~ de 

JV-1 

= ^2 (D2h(t k , q k ,t k +i,qk+i)Sq k + DiL(tk, q k ,t k+ i, q k+ i)5q k+1 )(t k+1 - t k ) 

fc=0 

JV-1 JV 

= 22 D 2 h(t k ,q k ,t k+1 , q k+1 )(t k+1 — t k )5q k + DiL(tk-i,qh-l,tk) qk)(fk — tk-i)Sq k 
fc=o fc=i 

JV-1 

= {D 2 'L{t k ,q k ,t k+ ±,q k+1 )(t k+1 - t k ) + D 4 IL(tfc-i, g^-i, k-, <7fc)(*fc - t k -i))Sqk 
fc=l 

+ D 2 L{t , go,ii, qi)(tx - to)Sqo + Z>4L(tjv-i, gjv-i, ijv, gjv)(ijv — t N -i)Sq N . 



(2.5) 



We can see that the last two terms in (2.5) come from the boundary variations. Based 
on the boundary variations, we can define two one forms on M x M 



and 



91 (t k ,qk,tk+i,qk+i) = D 2 'L(t k ,q k ,t k+ i,q k+1 )(t k+1 -t k )dq k , (2.6) 



^h + (^k, Ik, tk+i, <?fc+i) — D4h(t k ,q kl t k+ i,q k+ i)(t k+ i ~t k )dq k+1 . (2.7) 



Here we have used the notation in p2| . We regard the pair (9^~,9^ + ) as being the 
discrete version of the Lagrange one form V L in (2.4). 

Now we parameterize the solutions of the discrete variational principle by the initial 
condition (<7o;9i), and restrict § to that solution space. Then Eq. (2.5) becomes 

dS(q , qi, ■■■ , 9jv) • (6qo, 5qi, ■■■ , Sq N ) 
= 9^(t ,qo,h,qi) ■ (Sq ,8qi) + 6»+ (tjv-i, Qn-i, tN, qN) ■ (5q N -i,Sq N ) 

= ^l(*o,?o,*i,3i) • (&?0)<fyi) + (® N ~ 1 )*oi{to,qo,ti,<ii) ■ (8qo,5qi)- 



(2.8) 
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From ( |2.8| ), we can obtain 



dS = ei~ 



The Eq. (gj^ holds for arbitrary N > 1. Taking N=2 leads to 



By taking exterior differentiation of (2.10), we obtain 
From the definition of 0£" and 9^ , we know that 



(2.9) 



(2.10) 



(2.11) 



y L +ei+ = d((t k+1 -t k )h), 



which means d9^ + 



-ddt- 



Defining 



cjI =: dfff = -de^, 



(2.12) 



we show from (2.11) that the discrete flow <3? preserves the discrete Lagrange two form 

Let us consider an example. For the classical Lagrangian 

L(t,q, q) = ^-V(q), 
we choose the discrete Lagrangian h(t k , qk, tk+i, qk+i) as 

2 



it n + \ 1 I qk + 1 ~ qk 

^ \ ife+i — tk 



V 



q k +i + qk 



The discrete Euler-Lagrange equation (2.3) becomes 

qk+i - qk qk - qu-x \ . V'(q k ){t k+1 - t k ) + V'(q k -i)(t k - t k -i) 



. tk+i — tk tk — t k ^i 
which preserves the Lagrange two form 



= 0, 



1 



t k +i — t k 



fc+i k V"(q k ) I dq k+1 A dq k , 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



where q k = qk+q 2 h+1 . 



- _ qk-i+qk 

qk-i - — 2 — • 



If we take fixed variables tk+i — tk = tk — ifc-i = h, then (2.16) becomes 
q k+1 - 2q k + q k -i V'(q k ) + V'(q k -l ) 



h 2 



which preserves the Lagrange two form 



= 0. 
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2.2 Discrete mechanics based on finite element methods 

Now we consider discrete mechanics based on finite element methods. Let us go back 
to the variation problem of the action factional ([hi]) ■ The finite element method is an 
approximate method for solving the variation problem. Instead of solving the varia- 
tion problem in the space C 2 ([a,b\), the finite clement method solves the problem in a 
subspace Vftm([a, b]) of C 2 ([a, b]). Vh^([a,b]) consists of piecewise m-degree polynomials 
interpolating the curves q(t) S C 2 ([a,b]). 

First, we consider the piecewise linear interpolation. Given a partition of [a, b] 

a = t < t\ < ■ ■ ■ < tk < ■ ■ ■ < t/v-i < ijv = b, 

the intervals Ik = [*&, tfe+i] are called elements, hk = ifc+i — V/j([a,6]) consists of 
piecewise linear function interpolating q(t) at (tk, qk), k = 0, 1, • • • , N. Now we derive 
the expressions of qh(t) € Vh([a,b]). First we construct the basis functions fk{t), which 
are piecewise linear function on [a, b] satisfying (pk(U) = 51, i, k = 0, 1, • ■ • , N. 

0, otherwise; 0, otherwise; 



and for k = 1,2, • • • ,N - 1, 



1 + i fe -i^<^< fe ; 



Wfe(*) = ^ 1 _ t fc <t<tfc+i; (2-19) 

0, otherwise. 

Using these basis functions, we obtain the expression qh G Vh([a,b]): 

N 

q h (t) = J2lkfk{t). (2.20) 

fc=0 



In the space V/,([a, &]), the action functional (1.1) becomes 
S((t,q h {t))) = { L{t,q h (t),q h (t))dt 



N - 1 rt h+1 ( N N \ 

k=0 ^ tfc V i=0 i=0 / 

N-1 

^ L(i fc ,g fc ,t fc+ i,(7 fe+ i)(t fc +i -ife), ( 2 -21) 



fc=0 



where 

L(*fe, qk,tk+i,qk+i) 



tk+i — 



1 /•**+* / ^ d " \ 

— — / Lit, 22(qi(pi(t)), -S2(q^(t)) dt 
1 k •'tk \ i=0 j = o / 

ftk+l ( fc +! d fc +! \ 

y I ^(w(t)), s ^feViW) dt. (2.22) 



-k i—k 

Therefore, restricting to the subspace Vfc([a, b]) of C 2 (\a , b\), the original variational 
problem reduces to the extremum problem of the function (2.21) in qk, k = 0, 1, • • • , N . 
Notice that ( [2.21 ) is just one of the discrete actions (2.1). Thus, what remains to be 



G 



done is just to perform the same calc ulat ions on ( 2.21 ) as on (2.1). We can then obtain 
the d iscrete Euler-Lagrange equation (2.3) that preserves the discrete Lagrange two form 
( 2.12 ). Therefore, discrete mechanics based on finite element methods consists of two 
steps: first, use finite element methods to obtain a kind of discrete Lagrangian; second, 
use the method of Veselov mechanics to obtain the variational integrators. 

Let us consider the previous example again. For the classical Lagrangian (2.14), we 
choose the discrete Lagrangian li(tk,qk,tk+i,Qk+i) as 



1 




N 



\i=0 



tk+1 — tk 



1 ( Qk+i — qu 



V 



( tk+1 — t 



-qk 




\tk+l ~~ tk tk+1 — tk 



qk+i 



2 V tk+1 — tk 



- F{q k ,q k +i), 



where 



F(q k ,q^ 



1 



t k 



+i 



tk 



V 



tk-i 



t 



tk+i — tk 



qk 



t-t k 



tk-i 



tk 



qk+i dt 



dt 



(2.23) 



(2.24) 



The discrete Euler-Lagrange equation (2.3) becomes 
'qu+i - qk q k - qk-i 



tk+1 ~ 



tk — tk- 



dF( gk , qk+ i \ (th+i ^ + dF( q k-„ qk ) {tk _ ^ = Q) (2 ^ 5) 



dq k 



which preserves the Lagrange two form 



1 



tk+i ^ tk 



+ {t k+1 - t k ) 



d 2 F(q k ,q k +l) 
dqkdq 



k+l 



dq k 



dq k +i A dq k - 



(2.26) 



Again, if we take fixed time steps t k +i — t k = t k — t k -i = h, ( 2.25 ) becomes 
q k +i - 2q fc + qk-i , dF(q k ,q k+1 ) dF(q k _i,q k ) 



h 2 dq k 
which preserves the Lagrange two form 



1 h d 2 F(q k ,q k+1 ) 
h dq k dq k +i 



dqk 



dq k+1 A dq k - 



= 0, 



Suppose q k is the solution of ( 2.2q ) and q(t) is the solution of 

fq dV(g) _ 
dt 2 + 8q 

then from the convergence theory of finite element methods || [IlJ , we have 

\\q(t)-q h (t)\\^Ch 2 , 

where || • || is the L 2 norm. qu{t) = SfcLo Q k > h = m ax{/ifc} and C is a constant indepen- 

k 

dent of h. 



(2.27) 



(2.28) 
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If we use mid-point numerical integration formula in (2.24), we obtain 



i r h+1 Tr ( t k+ i-t , t-t k 

- tk J tk 

Qk + qk+i 



F(q k ,q k+1 ) =- — / V — q k + — <?fc+i ) 

In this case, is j ust (|2~T^ ). We can also use trapezoid formula or Simpson formula 

and so on to integrate (2.24) numerically and obtain another kind of discrete Lagrangian. 

2.3 Discrete mechanics with Lagrangian of high order 

Now we consider piecewise quadratic polynomial interpolation, which will result in a kind 
of discrete Lagrangian of high order. To this aim, we add an auxiliary node t k+ i to each 
element I k — [tk, t k +i}. There are two kinds of quadratic basis functions: 4> k (t) for nodes 
t k and 4> k+ i{t) for t k+ i that satisfy 

<t> k (t i )=5 1 ?, ^ fc (t i+ i) = 0, 

<t> k+i (t i+ i) = 5*, </> k+ i(ti) = 0, »,* = 0,1,.-.,JV. 
We list the basis functions as follows: 

fcW-K^-OOsf 1 - 1 )- (2.29) 
{ 0, otherwise; 

(j) N (t) = l\ hN -' J\hN-i J' (2.30) 
{ 0, otherwise; 

and for k = 1,2, • • • , TV - 1, 

0, otherwise; 
and for k = 0,1,--- ,JV- 1, 

^ fc+ i(t)= V *^ fc+1 ' (2.32) 

[ 0, otherwise. 

Using these basis functions, we construct subspace Vj,2([a, b]) of C 2 ([a,b]): 

N N-1 

q h 2(t) =Y,qMt) + 1k+i<t>k+i(t), Qh»(t) G V h »([a,6]). (2.33) 

fe=0 k=0 

In the space Vj,2([a, 6]), the action functional ( fL~l| ) becomes 



S((t,q h i(t)))= / L(t,q h 2(t),q h 2(t))dt 

J a 

^-1 /•*«,+ ! 



£ / " +1 L(t,q h2 (t),q h 2(t))dt 

k=0 , ' tk 
N-1 

2J L(i fc ,g fe ,g fc+ i,t fe+ i,% + i)(f fe+1 - i fc ), (2.34) 
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where 



L(£fc, Qk, 1 , tk+i, qjt+i) 



tk + l 



L{t,q h 2 (t),q h 2 (t))dt. 



tk+l - h Jt k 
For the discrete action ( |2.34 ), we have 

dS(q ,qi ,gi, • • ■ , q N _ 1+ 1 , qn) ■ (Sq ,Sqi,6qi, ■ ■ • ,Sq N _ 1+ i,8q N ) 

JV-l 

= J2 (DzL(w k )5q k + D a L{w k )Sq k+ i + D 5 h(w k )Sq k+1 )(t k+1 - t k ) 

fc=0 

N-l N 

= ^2 (D2h(w k )5q k + D 3 'L(w k )Sq k+ i)(t k+ i - t k ) + } D$L{wk-i){tk - i fe _i)<5g fc 



(2.35) 



k=0 
JV-l 



^2 (D 2 lL(w k )(t k+1 - t k ) + D 5 TL(w k -i)(t k - t k -i))5q k + D 3 TL(w k )Sq k+ i(t k 



fc=i 



+ £> 2 L(«>o)(ii - to)5qo + D 5 h(w N -i)(tN - t N -i)SqN, 



-t k ) 
(2.36) 



where io/c = (t k ,q k ,q k+ i,tk+i,q k +i), k = 0,1, ••■ , N - 1. From ( |2.36| ), we obtain the 
discrete Euler-Lagrange equation 



D 2 Uw k )(t k+1 - t k ) + D 5 L(w k -x)(t k - tfe-i) = 0, 
D 3 h(t k ,q k ,q k+ i,t k+ll q k+ i) = 0, 

Z3 3 L(tfe_i,gfe_i,g fe _ 1+ i,f fc ,%) = 0. 



(2.37) 
(2.38) 
(2.39) 



We can fro m (|2.38j ) and (2.3E) solve for and q k _i + i respectively, then substi- 

tute them into ([2.37| ) and finally solve for q k +\. Therefore, the discrete Euler-Lagrange 
equation ( 2.37 - 2~39| ) still determines a discrete flow 

* : M x M -> M x M 
^/(t k -i,q k -i,t k ,q k ) = (t k ,q k ,t k+1 ,q k+ i). 

From ( |2^3l) , we can define two one forms 

Ol~(t k ,q k ,q k+ i,t k+ i,q k+ i) = D 2 L(t& , <? fc , q k + A , *fc+l > ?fc+l) (*fc+l ~ *fc) d »> 

and 

0]L + (*Ai)3Ai,gjfe+i,tfc+l,gfc+l) = £ , 5M*fc,gfc,Q'fe+i,tfc+l)Q'fc+l)(tA!+l - t k )dq k+1 . 
Using the same method as before, we can prove that 

**{dQl + ) = -del'. (2.40) 



From the definition of 8^ and , we have 



L + L + = d ((tk+l - t k )L) - D 3 L(t k ,q k ,q k+ i,t k+1 ,q k+1 )dq k+ i. 



(2.41) 



From ( |2.38| ), we obtain D 3 L(t k , q k , q k+ i,tk+i,q k +i) = 0. Thus, 

©I" + ©I + = - t*)L), 

which means 



der = -d&t 



(2.42) 
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From fl2.4C| ) and (|2.42[ ), we arrive at 



(2.43) 



where fi£ = d6£ + . 



For the classical Lagrangian (fuij ), From ( [Egf ) and ( ^35| , we obtain 



-(g, l2 (t)) 2 -y(^(t)) )dt 



\ Q« 2 (4+i + + 4) + ab (tk + t k+l ) + b 2 



(2.44) 



where 



a= J{2 [ik+Qk+i - 2 1k+±- j • 



b = ( 4 fe + ^+1)^+1 - ( 3i fc + £fc+i)<?fc+i - (ife + 3t k+1 )q k 



and 



G(<7fc,g fc+ i,gfc+i) 



1 



tk+l ~ tk 



tk 



V (q k fk{t) + qk+ifk+i(t) + dt, 



where 



fk(t) 



= ( 2(t-t k ) Mt-tk 1 



t — tk f t — tk 

hk. 



For the discrete Lagrangian (2.44), the discrete Euler-Lagrangian equation (2.37-2.39) 
becomes 



aiqk-i + a 2 qk + »3<?fc+i + a 4<? fc -± + a^fc+i - d\hk - d 2 h k -i = 0, 



8 , 

^2 + -2g fe +i 



9G(q k) q k+ i,qk+i) 



dq 



o, 



9G(g fe _i,g fe _ 1+ i,%) 
9fc-i + ?fe - 2g fe _ 1+ 1 — = 0, 



(2.45) 
(2.46) 

(2.47) 



where 



1 1 



7 / 1 



ai = 



«4 



di = 



3 /ifc_i ' 
8 1 



a 2 = T, 



3 \h k -i h k 



8 1 



3ft.fc_i' 5 3 h k ' 



dG(q k ,q k+ i,q k+ i) 
dqk 



»3 



1 1 

3 ft-fc ' 



d 2 = 



dG(q k -i,q k _ 1+ i,q k ) 
dqk 
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The solution of (2.45-2.47) preserves the Lagrange two form 



1 d 2 G(q k ,q k+ i,q k+1 ) \ 

777 hk 5 — 77 M \dq k A dq k +i, 

3hk dq k dq k+ i I 



(2.48) 



where 



M 



J6. , d 2 G(Q k ,q k+ l,q k+ l) 
3hk k dq, , i dq k 



lfj d 2 G(q k ,q k+ i,q k + 1 ) 

3h- + ftfe 



9? fr , 1 9<?fc 



32 l 

37TT ~ 



2 G(<?fc:9 fc+ l ,q k + i) 
fc+ 3 



If we take the fixed time steps hk-x = h k = h, then (2.45-2.47)) become 
Qk-i - &q k - i + 14g fe - 8<? fc+ i + g fe+ i 



- dih k - d 2 hk-i = 0, 



8 / 

777^ (flk + fffc+x - 2?fc+I 



8 / 

777^ + 9fc - 2g fe _ 1+ i 



dG(q k ,q k+ i,q k+1 



0, 



dG(g fc _i,g fc _ 1+ i,g fc ) 



= 0, 



(2.49) 
(2.50) 

(2.51) 



which preserve 



' 1 d 2 G{q k ,q k+ i,q k+ i) \ 

77T - ,l a fl M dq k Adq k+ i, 

3h dq k dq k+ i 



(2.52) 



where 



M 



9 2 G(9fc,9 fc + l,9fc+i) 



JJ_ + h 

3fcfc 9q k+ idq k 



9 G(q k ,q k+ 1 ,9fc+i) 



16 + ft 

3ft 9g fc+ i 9gfc 



32 _ ^ 9 2 G(g^g fc+ i ^fc + i) 

3/1 ~~ ft a? - ; 



Suppose q k is the solution of (2.45-2.47)) and q(t) is the solution of (2.27), then from 
the convergence theory of finite element methods |5[ [llj , we have 



\\q(t)-q h 2(t)\\^Ch 3 



(2.53) 



where 



N 



N-l 



fc=0 



/i = max{/ifc} and C is a constant independent of h. 

k 



3 Time steps as variables 

In §2, the time steps t k play the role of parameters. They are determined beforehand 
according to some requirements. In fact, we can also regard t k as variables and the 
variation of the discrete action with respect to t k yields the discrete energy conservation 
law. This fact was first observed by Lee [ ||, ^, [l0| . The symplecticity of the resulting 
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integrators was investigated in @, [|] . These results also apply to the discrete mechanics 
based on finite element methods. 

We regard t k as variables and calculate the variation of the discrete action (2.1) as 
follows 



• • , tw, <?a0 • {Sto, Sqo, • • • , 5tN,SqN) 
S(to + eSt , q + e^o, ■ • ■ , t N + edt N ,q N + eSq N ) 



e=0 



dS(t , q 

d_ 
~ de 
N-l 

= ^ [D 2 L{w k )(t k+1 - t k ) + D 4 h(w k -x)(t k - t k -!)} 6q k 

k=l 
N-l 



^2 [DiL(wk)(tk+i - t k ) + D 3 L(w k -i)(t k - t k -i) + L(tOfc_i) - L(w k )] St k 



k=l 



+ D 2 L(w )(t 1 - t )Sq + DiL(w N -i)(tff - t N -i)6qN 

+ [DiL(too)(*i - t ) - L{w )] St 

+ [D 3 L(w N -i)(t N - <at_i) + L(«jjv-i)] St N , 



(3.1) 



where w k = {t k ,q k ,t k+ i,q k+ i), k = 0, 1, • • • , N - 1. From (|3.l[), we see that the 
variation Sq k yields the discrete Euler-Lagrange equation 



D 2 Uw k ){t k+ i - t k ) + DMw k -i){tk - tk-i) = 
and the variation Stk yields the discrete energy evolution equation 

Dih(tOk)(tk+i - t k ) + D 3 L(w k -i)(tk - tk-i) + L(wk-i) - L(w k ) = 0, 



(3.2) 



(3.3) 



which is a discrete version of ( |l.9| ). For a conservative L, (3^3) becomes the discrete 
energy conservation law. 

From the boundary terms in (|3.lD, we can define two one- forms 



~(w k ) = (DiL(«jfc)(tfc+i - tk) - L(w k ))dtk + D 2 1 L(w k )(t k+1 - t k )dq k , 



and 



9^{w k ) = (D 3 JL(w k )(t k+1 - tk) + h(w k ))dt k+1 + D 4 h(w k )(t k+1 - t k )dq kA 



(3.4) 



(3.5) 



These two one forms are the discrete version of t he e xtended Lagrange one form ( 1.12 ). 

Unlike the continuous cas e, the solution of (3.2) does not satisfy in general. 

Therefore, we must solve ( |3.2|) and (3.5) simultaneously. Using the same method in §2, 
we can show that the coupled integrator 



D 2 h{w k )(t k +i - t k ) + D 4 h(w k -i)(tk - t k -i) = 0, 

DiL(t» fc )(i* + i - t k ) + D 3 L(w k -i)(t k - t k -i) + L(w k -i) - L(w k ) = 

preserves the extended Lagrang e two form wl = d6^ . 
For the discrete Lagrangian (2.23), (3.6) becomes 



(3.6) 



( 



lk + 1 — Ik Ik — gfe-1 

tk — tk-1 
2 

+ F(w k 



V tk+1 ~ tk 

1 ( qk+i - qk 

2 \ tk+l ~ tk 



dF(w k ) dF(w k -i) h 

-Ilk H 7, llk-1 = U, 



dq k 

dF(w k 
dt k 



dq k 

qk — Qk-i 
tk — tk-l 



, , dF(w k -i) , 

+ F(w h -i) H tt: -h k -i. 

dtk 



For the kind of high order discrete Lagrangian, we can obtain similar formulas. 
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4 Conclusions 



In this paper, we have presented the discrete mechanics based on finite element meth- 
ods. The finite element methods provide a natural way to obtain a kind of discrete 
Lagrangian. Using piecewise quadratic functions, we have constructed a high order 
discrete Lagrangian and derived the corresponding variational integrator (the discrete 
Euler-Lagrange equation). The symplecticity of the high order variational integrator has 
also been proved. 

The time steps can play the role of parameters or the role of variables. We can obtain 
discrete energy conservation law when the time steps are regarded as variables. 

Finite element methods are well-established in numerical mathematics. Many results 
of finite element methods can be applied to the discrete mechanics presented in this 
paper such as convergence theory, error estimates as well as the solving method of the 
variational integrators. 

Recently, it has been proved ra] that the simplectic structure preserves not only on 
the phase flow but also on the flows with respect to symplectic vector fields as long as 
certain cohomological condition is satisfied in both continuous and discrete cases. This 
should be able to extend to the cases in this paper. We will study this issue elsewhere. 
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